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In order to understand why carbon dioxide (CO2) and methane (CH 4 ) molecules interact differ¬ 
ently with surfaces, we investigate the Casimir-Polder energy of a linearly polarizable CO2 molecule 
and an isotropically polarizable CH 4 molecule in front of an atomically thin gold film and an amor¬ 
phous silica slab. We quantitatively analyze how the anisotropy in the polarizability of the molecule 
influences the van der Waals contribution to the binding energy of the molecule. 

PACS numbers: 34.20.Cf; 42.50.Lc; 71.15.Mb 


I. INTRODUCTION 

Underground geological structures, e.g. rocks, shales 
and depleted coal beds, may contain large fractions of 
the world’s future fossil energy in the form of light nat¬ 
ural gas like methane [Ij. New production methods like 
hydraulic fracturing or fracking by CO2 injection allow 
trapped hydrocarbons to be produced directly from ti^t 
source rocks such as shale gas and shale oil systems [3. 
The reason that this method works may be partly due to 
the difference in electron affinities of the CO2 and CH4 
molecules Q. Much basic research focuses on molec¬ 
ular physisorption/chemisorption and meso-scale trans¬ 
port processes in nanostructured shales Q- By under¬ 
standing the underlying physical processes, the ultimate 
aim of such research is to realize a controlled displace¬ 
ment of methane (CH4) by carbon dioxide (CO2) in hy- 
drofractured shale and other formations. The added ben¬ 
efit of such a process is the simultaneous sequestration of 
the CO2 gas. CO2 gas can be trapped from huge point 
sources like power plants before it gets released into the 
atmosphere, and can be conveniently processed and uti¬ 
lized for fracking purposes j;^. 

From previous studies ^ 01 which do not take 
anisotropy into account, we know that a CO2 molecule 
and a CH4 molecule adhere to surfaces with very similar 
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FIG. 1 . (Color online) (Left) Schematic figure of an anisotrop- 
ically polarizable molecule above a dielectric slab. (Right) 
Schematic figure showing anisotropy of CO2 molecule and the 
isotropy of CH4 molecule in their polarizabilities. 

van der Waals energies while observed geological phe¬ 
nomena indicate strong preference of CO2 molecule over 
CH4 molecule in surface adsorption. In the present work, 
we explore the physical aspect of the underlying mecha¬ 
nism in order to understand the difference in interaction 
behavior of CO2 and CH4 molecules with surfaces with 
the focus being anisotropy in the electric polarizabilities 
of the respective molecules (see FIG.[T]). We address this 
by studying the interaction energy of an anisotropically 
polarizable molecule in front of a dielectric slab of finite 
thickness which is anisotrmic only in the direction per¬ 
pendicular to the slab [^, 0 | . We study the contribution 
from anisotropic polarizabilities to the Casimir-Polder in¬ 
teraction energy of a CO2 and a CH4 molecule in front of 
atomically thin gold films and an amorphous silica slab. 
We choose amorphous silica in particular because sur¬ 
faces generated by hydraulic fracturing are mostly min¬ 
eral surfaces like amorphous silica. Amorphous silica is 
isotropic in nature. We explore the data of the dielec- 
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trie function of gold available in Ref. 0 to study the 
variation of interaction energy with film thickness incor¬ 
porating the effects of anisotropy of the gold film. 

In the section II, we present the formalism of the 
Casimir-Polder interaction energy between a completely 
anisotropic molecule and a dielectric slab which is 
anisotropic in the direction perpendicular to the surface. 
In section III, we briefly summarize the method used for 
the calculation of the dielectric properties of the slabs. 
We also briefly describe the procedure used to obtain the 
anisotropic polarizabilities of the molecules. The dielec¬ 
tric functions of amorphous silica and gold are based on 
the density functional theory, and the anisotropic polar¬ 
izabilities of CO 2 and CH 4 molecules are also obtained 
from ab initio calculations pj]. fl^. Together the dielec¬ 
tric properties thus obtained are used to determine how 
the difference in the nature of polarizabilities of CO 2 and 
CH 4 distinguish their interaction energies. We present 
our results in section IV, and end with a few conclusions 
in section V. 


Our configuration with an anisotropic molecule above 
a dielectric slab with isotropic polarizability in the x- 
y plane renders the interaction energy independent of 
(j). Thus, if necessary, one can choose (j) = 0, but we 
will not bother to do so. The Casimir-Polder interaction 
energy, in this case, at zero temperature in the Fourier 
transformed space is, 


E = —he 


00 277 

dC f k±dk± r d4>k e 
c ./ 27r ./ 27r 


— 2K,a 


2k 


-!{<), ( 5 ) 


which is a generalization of the result given in Ref. Q. 
The details of the derivation leading to Eq. m has been 
omitted for brevity. Here, 


I{i() = r^[K^{k± • a • ku) + • a • z)] 

-r-®C^[(z X ku) ■ a • (z X k_L)]- (6) 


II. FORMALISM 

Consider an anisotropically polarizable molecule de¬ 
scribed by a frequency dependent molecular polarizabil¬ 
ity 


a{uj) = ai(w)eiei -|- a2(a;)e2e2 -I- a3(w)e3e3, (1) 

at a distance a above an anisotropically polarizable di¬ 
electric slab of thickness d described by dielectric permit¬ 
tivity 

£{u}) = ( 2 ) 

where T-components are in the x-y plane containing the 
dielectric slab and ||-component is normal to the surface 
of the slab. See FIG. [T] For an isotropic material like 
amorphous silica, we can set £-*- = sH. Magnetic perme¬ 
abilities for both the molecule and the dielectric slab are 
set to 1 . 

Here, the principal axes of the molecule are 


ei = cos 136 + sin Pp), 

(3a) 

02 = — sin /3 0 - 1 - cos (3 p, 

(3b) 

63 = f, 

(3c) 


where (3 is the rotation about the unit vector 63 , and 
f, 6, and (f> are the unit vectors in the spherical polar 
coordinates, 

f = sin 9 cos (j)'k + sin 9 sin (py + cos 9 z, (4a) 

6 = cos 9 cos + cos 9 sin (py — sin 9 z, (4b) 
0 = — sin^x-I-cos^y. (4c) 


The particular choice of (kj^, z x kj_, z) basis facilitates 
separation of TM and TE modes. Specifically, 


kj_ = cos (jifcX -I- sin (pky, 
z X kx = — sin (()fex -|- cos pty- 


(7a) 

(7b) 


and are the reflection coefficients for TM- and 
TE-modes: 


r 


H 


r 


E 






(1-6-2'=®'^) 

1 _ 7 K^-K, \ -2K^d 

\k^+k. j ^ 


( 8 a) 


( 8 b) 


where 
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and, d is the thickness of the dielectric slab. The inter¬ 
action energy after performing the pk integration is 
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A. Perfect conductor limit 


In the perfect conductor static limit, when —> 1 

and —>■ — 1 and the molecular polarizability gets a 

contribution only from the zero frequency, we reproduce 
the known result for the Casimir-Polder energy between 
an anisotropic molecular and a perfectly conducting slab: 


E = —he 
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where last term in Eq. m uniformly integrates to zero. 
Thus in this case the orientation of the molecule has no 
effect on the interaction energy. 


B. Temperature dependence 


FIG. 2. (Color online) The perpendicular and parallel dielec¬ 
tric constants for N=l, 3, 6 and 15 atomic layers of gold, and 
for an amorphous silica slab (written as a-Si 02 in the figure) 
in terms of the Matsubara frequencies. The perpendicular 
components of the 3, 6 and 15 atomic-layer thick gold films 
almost overlap. The dielectric constants at zero frequency are 
shown on the //-axis. 1 au=6.57968 x 10^® Hz=27.212 eV. 


To account for the temperature (T) dependence, we 
simply replace the integration over imaginary frequen¬ 
cies^^ a summation over discrete Matsubara frequencies 
Cn Em, 


poo „ —2«:a 

E = -2kBTy^' dkBk^^—I{lCn). (12) 

n=0 -^0 


where (n = 2'KkBTn/h, kB is the Boltzmann constant, 
and the prime indicates that the n = 0 term should be 
divided by 2. I{iCn) is given by Eq. ([H]), with ( replaced 
by Cn- 


C. Non-retarded limit 


where we have suppressed the frequency dependence. 
The orientation dependence appears only in the last 
term, which vanishes for 0 = 0,7r. 

I 


In the non-retarded limit = 0, the finite tempera¬ 
ture interaction energy between an anisotropic molecule 
at a distance a above an isotropic half-space Eq. m 
turns out to be. 
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using A = Iq^j- and 0. This energy is pro¬ 


portional to 


III. DIELECTRIC FUNCTION AND 
POLARIZABILITY 


All calculations for amorphous silica were carried out 
using the Vienna ab initio simulation package (VASP) 
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FIG. 3. (Color online) The anisotropic polarizabilities of CO 2 
and CH 4 in units of A® in terms of the Matsubara frequencies. 
The zero frequency polarizabilities are indicated on the j/-axis. 
Note that CO 2 is much more polarizable in the transverse 
direction. 

with the Perdew-Burke-Ernzerhof (PBE) [T^ functional. 
Projector augmented wave (PAW) pseudopotentials [l^, 
[13 were used to model the effect of core electrons. The 
non-local parts of the pseudopotentials were treated in 
the real space for the Born-Oppenheimer molecular dy¬ 
namics (BOMB) and in the reciprocal space for all other 
density functional theory (BET) calculations. The struc¬ 
ture of amorphous silica was generated using the BOMB 
simulations of 72-atom supercell with different annealing- 
quenchjM temperature protocols similar to earlier stud¬ 
ies (l§ . [ 1 ^ . The dielectric properties of amorphous silica 
were then calculated using the scissors-operator approx¬ 
imation (A=3.6) for PBE calculations. The dielectric 
function on the imaginary frequency axis was determined 
using the Kramers-Kronig dispersion relation. The low- 
energy spectra are verified by calculating the static di¬ 
electric constants from the Born effective charges. The 
static dielectric constant was found to be 4.08±0.11. The 
details of the calculations of anisotropic dielectric func¬ 
tions for gold sheets were presented by Bostrom et al. in 
Ref. [ 13 . We plot the parallel and perpendicular dielec¬ 
tric constants as defined in Eq. (2) of different thicknesses 
of gold films, and of amorphous silica in FIG. [5] 

The anisotropic polarizability tensors at imaginary fre¬ 
quencies for CO 2 and CH 4 were calculated using the 
quantum chemistry package MOLPRO [13 . Calculations 
were performed at the coupled clusters, singles and dou¬ 
bles (CCSB) level of theory. The correlation-consistent 
aug-cc-pVQZ basis set (13 was used. The geometries of 
the molecules were first optimized by energy minimiza¬ 
tion before being used in polarizability calculations. All 
calculations were done at room temperature. We show 
the anisotropic polarizabilities of CO 2 and CH 4 molecules 

in FIG. [3 


FIG. 4. (Golor online) Comparing the interaction energy of 
CO 2 and CH 4 molecules at a distance of 10 A from an amor¬ 
phous silica slab with respect to 0-orientation. We refer to 
the configuration 0 = 0 as the parallel orientation and 0=7r/2 
as the perpendicnlar orientation. All energies are in units of 
keT. 



FIG. 5. (Color online) Interaction energy curves of a CO 2 
molecule for different 0-orientations near gold films of N- 
atomic-layer thickness. The corresponding energies for a CH 4 
molecule are -57.9, -77.8, -84.8 and -87.8 in units of 10“® ksT 
for N=l, 3, 6 and 15 respectively. 


IV. NUMERICAL RESULTS 

For a linear molecule like CO 2 , the two most notable 
configurations are the parallel and perpendicular orien¬ 
tations with respect to the dielectric slab, and we are 
interested in the change in interaction energies in going 
from one orientation to the other. By parallel orienta¬ 
tion, we refer to the configuration in which 63 is aligned 
along z while perpendicular orientation refers to the case 
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FIG. 6 . (Color online) Comparing the interaction energy of 
CH 4 and CO 2 molecules at the parallel {9 = 0) and the per¬ 
pendicular {6=n/2) orientations at a varying distance from an 
amorphous silica slab. The interaction energy for a methane 
molecule is independent of ^-orientation. The inset figure 
shows the small distance limit. The axis labels and the leg¬ 
ends of the outer figure hold for the inset figure as well. 



FIG. 7. (Color online) Comparing the interaction energy of 
CH4 and CO 2 molecules at the parallel {9 = 0 ) and the per¬ 
pendicular {9=tt/2) orientations at a varying distance from a 
15-atomic-layer thick gold film. 

when they make an angle 0 = 7 r/ 2 . A CO 2 molecule has 
anisotropy in one direction in its diagonal basis as shown 
in FIG. [ 3 l Thus, choosing 03 along the unique linear axis 
of the molecule, it is obvious from Eq. OT that there is 
no /3 dependent term in the interaction energy. In this 
particular choice of axes, the unique linear axis of the 
CO 2 molecule is perpendicular to the surface when 9 = 0 
(parallel orientation) and parallel to the surface when 
6=7r/2 (perpendicular orientation). The curves in FIG.|3] 


show the Gasimir-Polder interaction energies of a CO 2 
and a GH 4 molecule for different 6 *-orientations placed at 
a distance of 10 A from an amorphous silica slab. As ex¬ 
pected, a methane molecule being highly isotropic shows 
no change in energy with change in 9. A CO 2 molecule, 
on the other hand, exhibits a slight change in the in¬ 
teraction energy at different orientations. The curve in 
FIG. m shows that the GO 2 molecule has lower energy at 
the parallel orientation {9 = 0) than at the perpendic¬ 
ular orientation (0=7r/2) near an amorphous silica slab. 
Thus, the molecule is most stable when its unique linear 
axis is aligned perpendicular to the slab. Irrespective of 
their orientations, the magnitude of the energy is larger 
for a GO 2 molecule than for a GH 4 molecule. This may 
play a role in the observed preferential adsorption of CO 2 
over GH 4 molecules on surfaces. These interaction ener¬ 
gies are, however, very small compared to They 

become comparable when the molecule is very near the 
surface (see insets of FIGs. |n]and[71). It should be noted 
that the dielectric continuum picture of our model breaks 
down at the small distance limit (roughly below lOOA). 
A more rigorous quantum chemistry calculation will be 
required to take into account the effects due to surface, 
bonding, etc. 

To make an estimate within our model, we calculate 
the energy of a system consisting of an amorphous sil¬ 
ica slab with CO 2 molecules in the parallel orientation 
at, say, 8 A and CH 4 molecules at, say, 5 A from the 
slab mimicking the condition when the CO 2 gas is be¬ 
ing injected. We then consider the reverse system when 
the CO 2 molecules are at 5 A and CH 4 molecules at 8 
A. The difference in the interaction energies between the 
two configurations is Q.078kBT, which is roughly 18 % 
compared to the energy in the first configuration. Thus, 
the second system with CO 2 near the surface is more 
favorable. As stated earlier, at such small separation dis¬ 
tances there would be considerable contributions to the 
interaction energy from other effects. 

In FIG.O we plot curves for the variation of molecule- 
surface interaction energy with respect to 9 for a GO 2 
molecule near gold films of different thicknesses. As can 
be observed from the figure, thicker films give larger mag¬ 
nitudes of interaction energies. Only the interaction en¬ 
ergy with the 1 -atomic-layer thick gold film displays ap¬ 
preciable difference in comparison with the energy curve 
for N=15 atomic-layer thick gold film while the interac¬ 
tion energies with the 3 and 6 -atomic-layer thick gold 
films gradually approach that of 15-atomic-layer thick 
gold film. The energy corresponding to N=15 atomic 
layers of gold is largest in magnitude. From FIGs. 0] 
and [5l we can see that the trends of the energy curves 
are alike but the molecules have energies larger in mag¬ 
nitude for the more dielectric gold film compared to 100 
A thick amorphous silica. We also provide the corre¬ 
sponding energies for a CH 4 molecule near gold films of 
varying thickness in the caption of FIG. 01 They are 9 
independent. 

In FIG. 01 we fix the CO 2 molecule in the parallel and 
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perpendicular orientations and plot the interaction en¬ 
ergy with respect to separation distance from the amor¬ 
phous silica slab on a logarithmic scale. The energy curve 
for CH 4 , which is orientation-independent, is also shown. 
A small difference in the energy curves for the paral¬ 
lel and perpendicular orientations is observed for a CO 2 
molecule. The interaction energy is larger in magnitude 
for a CO 2 molecule than for a CH 4 molecule at all sep¬ 
aration distances from the slab owing to greater polar¬ 
izabilities of CO 2 molecule. The curves follow the 
dependence of the non-retarded approximation up to a 
separation distance of few Angstoms, and gradually de¬ 
viates. The inset figure shows the interaction energy in 
the small molecule-slab separation distance limit (on a 
linear scale). FIG. [7] shows similar curves near a 15- 
atomic-layer thick gold film. 

V. CONCLUSIONS 

In this work, we presented a generalized expression for 
the interaction energy between a completely anisotropic 
molecule and a dielectric slab polarizable in the direction 
perpendicular to the surface. Applying this to the spe¬ 
cific case of a linearly polarizable CO 2 molecule and an 
isotropically polarizable CH 4 molecule, we showed that 
the anisotropy influences the van der Waals energy to a 
small degree. The parallel orientation (0 = 0) is more 
favored in comparison to the (0=7r/2) perpendicular ori¬ 
entation in the case of a CO 2 molecule. In subsequent 


studies, it will be interesting to incorporate the effects 
of finite-size of the molecule in which one has to care¬ 
fully consider different radii of the anisotropic molecule 
in different directions for determination of the interaction 
energy for different orientations. In the future, we hope 
that it will prove possible to transcend the limitations of 
the continuum approximation, to get more reliable esti¬ 
mates of Casimir-Polder energies at very short distances 
than we can provide here. 
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